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Phase  Plane  Analysis  of  Dynamically 
Leaded  Journal  Bearings 


H.  T.  ALBACHTEN 


NOMENCUTURE 
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radial  clearance.  In. 
eccentricity.  In. 
nonlinear  perturbations 
a  constant,  l^Tr/^ir^L^ 

Wc2 

(  J  =  12  TT^S ) 
axial  length  of  the 
fluid  pressure  psl 
journal  radius.  In. 


Journal,  In. 


{x2 


-2,1/2 
+  y  ) 


=  bearing  radius.  In. 

=  Sommerfeld  number,  L  i 

JtWc2 


t  =  time,  sec 

W  =  magnitude  of  the  load  on  the  Journal, 
lbs. 

X  =  phase  plane  co-ordinate  =  t 

d  t 

y  =  phase  plane  co-ordinate  =  J  — 

dt 

x,y  =  co-ordinates  to  a  point  In  phase  plane 
measured  from  singular  point 
x*,y«  =  phase  plane  co-ordinates  of  a  singular 
point 

X(x,y)  =  denominator  of  right-hand  side  of  dif¬ 
ferential  equation  for  Integral  curves 
Y(x,y)  =  numerator  of  right-hand  side  of  differ¬ 
ential  equation  for  integral  curves 
z  =  x*^ 

fluid  viscosity,  lb  sec/ln.^ 

0  =  attitude  angle,  rad 

=  rotational  speed  of  Journal,  rad/sec 
*^2  =  rotational  speed  of  the  load,  rad/sec 
n  =  -  2  rad/sec 

f  =  eccentricity  ratio,  e/c 
o  =  subscript  denoting  Initial  value 
•  =  differentiation  with  respect  to  time 


motion  that  included  certain  previously  neglected 
second-order  terms.  In  1949  Shaw  and  Macks  (5) 
gave  a  rather  complete  up-to-date  summary  of  the 
subject,  as  well  as  an  extensive  bibliography. 

The  purpose  of  this  paper  Is  to  analyze 
Harrison's  equations  by  using  phase-plane  methods 
(resulting  In  velocity  versus  displacement  plots 
of  the  Journal  eccentricity)  and  by  this  treat¬ 
ment  acquire  greater  Insight  into  Journal-bearing 
behavior.  It  Is  believed  that  this  type  of  anal¬ 
ysis  exhibits  previously  known  Journal-bearing 
characteristics  more  clearly  and  Introauces  an 
explanation  of  certain  other  characteristics 
which  were  previously  obscure. 

For  example,  all  possible  motions  of  the 
Journal  that  can  result  from  varying  such  param¬ 
eters  as  Journal  rotational  speed,  load,  eccen¬ 
tricity,  and  so  on,  are  clearly  shown.  Thus,  the 
sensitivity  of  the  Journal  to  the  half-frequency 
Instability  Is  btought  out  by  the  insight  gained 
in  studying  the  effects  on  the  phase-plane  plots 
of  varying  these  parameters. 


Fig.  1  Journal  bearing  configuration  relative 
to  fixed  point  in  space 


2 

In  1920,  Harrison  made  the  first  dynamic 

analysis  of  the  Journal  under  the  Influence  of 
external  load  and  fluid  pressure  forces.  His 
classical  work  results  In  two  well-known  equa¬ 
tions  of  motion  ^(refer  to  equations  (1)  and  (2)^ 
from  which  Important  fundamental  Information  has 
been  derived.  Swift  (^)  and  Burwell  (2)  applied 
Harrison's  results  to  several  Important  examples, 
and  In  1933  Robertson  (4)  derived  equations  of 

2  Underlined  numbers  In  parentheses  designate 
References  at  the  end  of  the  paper. 


DESCRIPTION  OF  GENERAL  PROBLEM 
Equations  of  Motion 

The  system  considered  in  this  paper  Is  Il¬ 
lustrated  In  Flg.l.  The  dynamic  equivalent  of 
this  system,  shown  In  Pig. 2,  regards  the  load  as 
fixed  and  both  the  bearing  and  Journal  as  rota¬ 
ting. 

The  perfectly  aligned.  Infinitely  long,  mass¬ 
less  Journal  rotating  In  a  complete  bearing  while 
under  the  Influence  of  fluid  pressure  and  exter¬ 
nal  load  forces  was  analyzed  by  Harrison  who  de- 


rived  the  equations  (reference  5,  p.  226) 

T  d  I 

TTl/z 


cot  i  s 


dt 


(  I 


■  in  4 


J«  (Q  -2  II-) 


(1) 


(2) 


(  Z*  *h  n/i  -«  ^ 

where  the  eccentricity  ratio  t  =  e/c,  fS  is  the 
attitude  angle  measured  between  the  load  and  ec¬ 
centricity  vectors,  Tl  =  u-.^  -  2  ^'2  and  J  is  the 
constant 


J  = 


12  Ml 


(3) 


W 


J  =  12  TT^S  Where  S  is  the  Sommerfeld  number. 
When  0  varies  with  time,  d0/<it  is  the  whirl  speed 
(see  Nomenclature  for  the  definition  of  other 
terms ) . 

If  equations  (1)  and  (2)  could  be  solved  si¬ 
multaneously,  the  Journal  motion  would  be  deter¬ 
mined  by  the  solution  functions  ^(t)  and  t  (t). 
Since  this  does  not  seem  possible,  two  alterna¬ 
tive  procedures  would  be  to  either  consider  the 
four-dimensional  phase  space  in  the  displacements 
and  velocities  of  the  two  variables,  or  to  con¬ 
sider  simultaneously  two  separate  phase  planes. 

As  the  first  step  towards  the  latter  approach, 
equations  (1)  and  (2)  may  be  combined  and  0  elim- 


Fig.  2  Journal  bearing  configuration  relative 
to  fixed  point  on  the  load  vector 


inated  to  make  possible  the  derivation  of  the 
phase  plane  equations  for  the  6  variable  (see 
Appendix  1),  but  It  does  not  seem  possible  to 
eliminate  t  between  equations  (1)  and  (2)  In  or¬ 
der  to  form  the  phase-plane  equations  for  the  0 
variable.  However  the  solution  0  =  0{t)  may  be 
obtained  by  eliminating  dt  between  equations  (1) 
and  (2)  (reference  5,  p.  358),  which  would  make 
possible  a  polnt-by-polnt  construction  of  the  0 
phase  plane  from  the  (  phase  plane. 

In  general.  If  the  €  phase-plane  nonordinary 
points  are  first  determined.  It  Is  possible  to 


gain  Important  information  about  the  complete 
system  by  e.xaminlng  the  behavior  of  the  deriva¬ 
tives  d0/dt  and  d^0/dt^  at  those  points.  As  will 
be  seen  later,  this  general  approach  is  simple  to 
carry  out  for  this  particular  system. 

Singular.  Critical,  and  Equilibrium  Points 

Since  the  procedure  followed  here  will  be  to 
describe  the  Journal  motion  from  concurrent  in¬ 
vestigations  of  the  6  phase  plane  and  equations 
( 1)  and  (2),  it  is  essential  to  distinguish  be¬ 
tween  a  mathematical  singular  point,  a  critical 
point,  and  a  Journal  equilibrium  point.  This  is 
because  neither  a  singular  point  nor  a  critical 
point  for  the  one-degree-of -freedom  i  phase  plane 
is  necessarily  an  equilibrium  point  for  the  two- 
degi-ee-of -freedom  Journal  motion.  For  other  phe¬ 
nomenon  completely  descrlbable  by  one-degree-of- 
freedom,  critical  points  are  equilibrium  points. 

Possible  Indefiniteness  can  also  occur  if  it 
is  not  realized  that  although  insight  useful  for 
sketching  the  phase  plane  "Integral  curves"  (geo¬ 
metrical  configurations)  may  be  achieved  by  in¬ 
vestigating  the  mathematical  singular  points  of 
the  phase-plane  equation  written  as 

^ 

X  (  x,  y) 

some  of  these  points  may  not  even  be  critical 
points  of  the  system  written  as 

X  =  ^  =  y  Y  =  4^  =  $(  X,  y  )  ^ 

Which  define  the  phase-plane  "trajectories."  In 
this  terminology  the  Integral  curves  of  equation 

(4)  are  the  same  as  the  trajectories  of  equation 

(5)  but  lack  their  "time  sense"  (reference  6,  p. 
116) .  The  term  "singular  point"  will  be  used  to 
designate  a  nonordinary  point  on  an  Integral 
curve  while  "critical  point"  will  refer  to  such 
a  point  on  a  trajectory.  Each  critical  point  of 
equation  (5)  is  usually  a  singular  point  of  equa¬ 
tion  (4).  In  all  cases  the  Journal's  equilibrium 
points  are  only  at  the  locations  where 

d«  /dt  =  d^«/dt^  =  d^J/dt  =  d^^  /  d  =  3 
If  the  following  derivative  expressions  are 
formed  from  equations  (1)  and  (2),  the  exact  lo¬ 
cation  of  the  equilibrium  points  will  be  made 

=  -j-  (  1  -  f  3/2  ^,^5  ^ 

=  ^  -  sin  i  (  2  -I-  «  ^)  \/i  -  €  ^ 

^  2  J« 

-  -J  •>" '  U 

,  J  M  1  V  ,  ii-j 


clear; 

d  t 
d  t 
d  i 
d  t 


2 


i  d  ^ 

cos  t  :5 — 

d  t 


_L_  [(  2  +  ‘  ■  ‘  ■ 

2J  [ - ^  - 

/2e^-f‘'t2  ,  ,,d«  1 


(6) 


Equations  (6)  show  that  only  the  point  located 
at  =  ti/2,  and  the  particular  value  of  i  which 
makes  d0/dt  =  o  is  an  equilibrium  point  for  the 

p 

system.  At  6  =  1  all  derivatives  are  zero  ex¬ 

cept  dfi/dt,  which  is  then  equal  to  ■^/2.  This  is 
the  half-frequency  whirl  phenomena  which  will  be 
the  instability  encountered  at  that  phase  plane 
critical  point. 

Phase-Plane  Equations 

The  differential  equation  defining  the  radial 
motion  of  the  Journal  center  is 


d 


(7*^+2)  ,d«  ,2 

2  .  (1  -  d  t 


(7) 


+ 


jn 

2J 


where,  in  accordance  with  the  notation  Introduced 
in  Appendix  1 ,  0-^0  ^  tt  and  £  is  positive  when 
the  journal  center  is  to  the  right  of  the  fixed 
vertical  reference  in  Pig. 2,  and  negative  when 
the  journal  center  is  to  the  left  of  the  vertical 
reference  line.  It  is  understood  that  this  nota¬ 
tion  is  a  mathematical  convenience  only,  since  a 
negative  £  has  no  physical  meaning. 

Equation  (7)  may  be  put  in  conventional,  non- 
dimensional  phase-plane  notation  and  form  by  let¬ 
ting  fc  =  X  where  -1  ^  x  ^  +1  and  J  d£/dt  =  y. 
The  equations  defining  the  €  phase-plane  trajec¬ 
tories  are  then 


d  X 

d  t 


d  y  _  -  (7  Z) 


d  t 


X^) 


o 


2  X  (1 

\/(i-xV-y== 


f 


(8) 


+  (1 


-x¥ 


(2  t 


2  X  J 


and  the  equation  defining  the  £  phase-plane  inte¬ 
gral  curves  may  be  put  in  the  following  form: 


Fig.  3  Ph  blot  for  A  =  0 


(8).  The  singular  points  of  the  Integral  curves 
are  at  the  locations  which  simultaneously  make 
X(x,y)  and  Y(x,y)  =  0  in  equation  (9)-  The  co¬ 
ordinates  of  either  type  of  point  are  denoted  by 
{x*,y*) .  The  co-ordinates  measured  from  such  a 
point  are  distinguished  by  a  bar.  Thus, 

X  =  X  -  X*  (10) 

7  =  y  -  y* 

A  procedure  (outlined  below,  but  which  is  ex¬ 
plained  in  detail  by  Coddlngton  and  Levinson, 
reference  7)  Tor  determining  the  nature  of  a  cri- 


d  y  _  -(7  x^  t  2)  -  jQxil 


X^)^  -y^+  (1  -  x^)  ^  (2+  x^)  _ 


d  X 


g  Y(x  y) 


2  X  y  (1  -  x‘) 


X(x,  y) 


(9) 


\  All  Integral  curves  in  the  phase  plane  as  defined 

by  equation  (9)  will  be  symmetrical  about  the  x- 
axis,  but  the  x  in  the  term  containing  the  square 
root  destroys  the  symmetry  about  the  y-axls. 

Evaluation  of  Nonordinary  Points 

Critical  points  in  the  t  phase  plane  are  val¬ 
ues  of  {x,y)  that  make  x  =  y  =  o  in  equations 


tical  point  is  to  first  transform  the  phase-plane 
equations  into  those  whose  co-ordinates  are  mea¬ 
sured  from  the  critical  point  (x*,y*)  by  using 
equations  (10).  The  system  then  considered  is 
one  where  this  transformation  results  in  expres¬ 
sions  of  the  form 


3 


y 


(a)  Phase  Plane  (b)  Journal  Motion 

Fig.  4  Configurations  for^l  =  o 


^  =  cx+dy+l,  Ot,  y) 

(11) 

=  ax  +  by+lj  (X,  y) 

Where  a.b.c.d  are  real  constants,  ad  -  be  ^ 0  and 
fj,  are  real  continuous  functions  defined 
within  some  circle  about  (x*,y*)  with  radius  r^ 

>  0.  The  functions  fj^  and  fg  are  called  pertur¬ 
bations  and  the  nonlinear  system  described  by 
equation  (11)  Is  referred  to  as  the  system  per¬ 
turbed  from  the  linear  system  given  by 

lx  =  c  X +  d y 

dt  (12) 

lx  =  a  X  +  b  y 

d  t 

After  the  conversion  of  f^  and  ^2  Into  polar  co¬ 
ordinates,  If  fi/ro  and  f2/ro  approach  zero  as 
Tq  — >  0  +  the  perturbations  are  then  guaranteed 
to  tend  to  zero  faster  than  the  linear  terras  In 
equation  (11).  Under  these  conditions  the  behav¬ 
ior  of  the  trajectories  of  equation  (11)  near  the 
critical  point  are  very  similar  to  the  behavior 
of  the  trajectories  of  equation  (12).  For  the 
two  cases  considered,  the  theorems  given  In  {J) 
are  referenced  freely  In  describing  the  Journal 
bearing's  phase  plane  portraits. 

THE  CASE  OF  fl  =  0 

Phase  Plane  Portraits 

For  the  special  case  of  the  load  rotating  at 
one  half  the  Journal  speed  (fl  =  -2  u.'2  =  0)  or, 

equivalently,  the  case  of  a  nonrotating  Journal 
acted  upon  by  a  nonrotating  load,  equation  (7) 
reduces  to  the  following: 

X  7*^+2  /  d  «  y  (I  {Z*  (13) 

2«(1  -  dt  '  2«J^ 

For  this  case.  It  Is  more  meaningful  to  Investi¬ 
gate  the  equation  for  the  Integral  curves  which, 
from  equation  (13),  becomes 


d  y  ^  -  (7  X* 

d  X  ~ 


+  (1  -  x^)  ^  (2  t  x^)  ^  Y  (x,y)  (14) 


2  X  y  (1  -  x‘) 


X  (x,  y) 


Where  now  the  phase-plane  portrait  Is  symmetrical 
about  both  the  x  and  y-axes . 

The  singular  points  In  equation  (14)  are  at 
the  four  locations  (+  1,  0)  and  (0,  +  1).  Al¬ 
though,  for  completeness,  a  brief  discussion  of 
these  singular  points  will  be  given  In  the  next 
section,  their  evaluation  Is  not  necessary  In  or¬ 
der  to  obtain  the  phase-plane  portrait.  This  Is 
because  equation  (13)  may  be  Integrated  once  to 
yield  an  explicit  equation  for  the  Integral 
curves  and,  when  the  constant  of  Integration  Is 
zero  (xq  =0),  It  may  be  Integrated  again  to  give 
X  as  a  function  of  time.  The  result  of  the  first 
Integration  (see  Appendix  2)  Is 


r  2  II  2.31  „  2.  9/2 
[y.  -  (I  ■ «.  I J  I  . 


(1  - 


Prom  equation  (15)  It  Is  seen  that  the  Integral 
curves  meeting  the  y-axls  will  be  given  by 

y=  +  ,l.x2)5/2  (16) 

Equation  (l6)  defines  the  boundary  of  a  closed 
area.  Each  boundary  arc  connects  singular  points 
and  Is  called  a  separatrlx  arc.  It  Is  easily 
seen  that  all  physically  meaningful  Integral 
curves  must  lie  Inside  the  boundai^  defined  by 
equation  (l6),  for  If  y  >  (l-x^)^/^  ory<- 
(1-x^)^'^,  then  equation  (1)  would  Imply  cos  fS'> 

1  or  cos  <  -1  .  Moreover,  the  difference  be¬ 
tween  the  ordinate  of  an  Integral  curve  beginning 
Inside  the  separa trices  and  the  ordinate  of  a 
separatrlx  Itself  remains  finite  and  tends  to 
zero  only  at  the  singular  points  (+  1,  0).  The 
Integral  curves  are  tangent  to  the  separatrlxes 
at  (+  1,  0),  but  they  do  not  touch  them  at  any 
other  point.  Fig. 3  shows  Integral  curves  plotted 
from  equation  (15). 

Under  zero  Initial  conditions  the  solution 
of  equation  (l6)  Is 


4 


X  = 


t 


(17) 


which  gives  the  change  in  £  with  time  along  a 
separatrix  arc  (trajectory). 

Examination  of  equations  (1)  and  (2)  reveals 
how  0  varies  with  time  along  a  separatrix  arc. 
For  convenience  express  equation  (1)  and  (2)  in 


phase-plane  notation,  thus 

cos  i  =  _ I _ , _ 

(1  -  x2)3/2 

(18) 

sin  d  =  j  X  (  0  -  2ar ) 

(19) 

(2  +  x^  )\/l  -  x^ 

When  fS  =  0  and  tt  equation  (l8)  reduces  to  equa¬ 
tion  (l6).  When  ji  =  0  or  It  equation  (19)  re¬ 
quires  dfi/dt  =  0  for  n  =  0  and  0  <  x  <1.  The 
separatrix  arc  between  (0,  +1)  and  (+1,  0)  is  the 
Integral  curve  along  which  0=0  and  d^/dt  =  0 
while  the  path  between  (1,0)  and  0,  -1)  requires 
0  =  Ti  and  A0/dt  =  0.  Similar  statements  apply  to 
the  negative  half  plane. 

The  physical  interpretation  of  Journal  behav¬ 
ior  can  be  most  easily  explained  by  referring  to 
Pig. 4.  As  trajectory  A,B,C,D  in  Pig. 4(a)  is  tra¬ 
versed  in  the  phase  plane,  tne  Journal  travels 
vertically  down  from  wall  contact  at  =  tt  to 
wall  contact  at  0  =  0.  For  other  trajectories 
such  as  A,J,K,D,  the  Journal  motion  is  to  the 
right  of  the  vertical  reference  as  Indicated  in 
Fig. 4(b).  Paths  are  to  the  left  of  the  vertical 
reference  for  trajectories  such  as  E,M,N,H,  This 
type  of  Journal  behavior  has  been  previously  com¬ 
puted  and  sketched  by  Swift  (reference  5,  p.  222). 


THE  CASE  OF  A  7^  0 

Singular  Points  in  the  Phase  Plane 

In  general,  and  the  integral  curves  are 

defined  by 


d  y 

d  X 


-  (7x^  +  Z)  -  JOx(1  -  x“)  Vfl  -  x^)^  -  +  (1  -  x^)^(Z  +  x^) 

2  X  y  (  1  -  x^) 


Y(x.y) 
X(x,  y) 


(9) 


Phase-Plane  Nonordinary  Points 

Among  the  four  nonordinary  points  of  equation 
(14),  only  those  at  (+  1,  0)  are  critical  points 
and,  from  equation  (6),  these  points  are  equilib¬ 
rium  points  for  the  system.  The  theorems  of 
Coddlngton  and  Levinson  (2.)  indicate  unstable 
"node-like"  critical  points  at  (+1,  0)  and  "sad¬ 
dle  like"  singularities  at  (0,  +  1). 

The  case  of  A=  0  is,  by  Itself,  of  only 
slight  Interest,  but  it  does  serve  as  a  founda¬ 
tion  for  the  more  Involved  case  of  A  5^0.  In  ad¬ 
dition,  it  gives  an  indication  of  the  Journal  be¬ 
havior  as  it  slows  down  for  stopping  (if  “^2  “ 
or  as  the  load  rotational  speed  approaches  one 
half  the  Journal  rotational  speed. 


which  has  singular  points  at  (0,  +1)  and  at  the 
points  along  the  x-axls  where 
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(2  +  x*^)  J  =  0  ( 


20) 


Equation  (20)  yields  the  singular  points  at  (+  1, 
0),  and  a  fifth  singularity  whose  x*  location  is 
contained  in  the  solution  of  the  cubic  equation 


+  3z^  +  j2q2  ^ 


4=0 


(21) 


Where  z  =  x*2.  Consideration  of  equation  (21)  in 
the  root  locus  form 


0  ^  z  -  4 

1  +  — z- 


z  (  z  +  3) 


(22) 


shows  that  it  will  always  have  a  positive  z  root 


(a)  Phase  Plane  (b)  Journal  Motion 

Fig.  6  Configurations  for  large  values  of  J  -iT- 


located  between  z  =  0  and  z  =  4/j2I12.  por  val¬ 
ues  of  X*  Is  closely  approximated  by 

X*  =  2/JQ  for  ^  6,  (23) 

l/2 

where  only  the  positive  sign  resulting  from  Z  ' 
satisfies  equation  (20).^  The  values  of  x*  sat¬ 
isfying  equation  (20)  are  plotted  In  Fig. 5  which 
Is  the  familiar  Sommerfeld  eccentricity  curve  de¬ 
fining  the  equilibrium  position.  Thus  the  singu¬ 
lar  points  foril^O  Include  those  foril_=  0  plus 
a  fifth  point.  Before  Investigating  the  nature 
of  these  singular  points.  It  Is  desirable  to  ex¬ 
amine  equations  (l8)  and  (19)  under  conditions 
when  n.  ft  =  0 


C08  i  =  y 

(1  - 

(18) 

8ln  d  =  J  X  (Cl  -  2  3r) 

(19) 

(2  +  x^)  Vl  -  x^ 

0  and  V  equation  (l8)  reduces 

to  equation 

(l6),  which  from  equation  (9),  Is  also  an  Inte¬ 
gral  curve  of  the  6  phase  plane  as  In  the  case  of 
A.=  0.  The  right  side  of  equation  (19)  can  then 
only  be  zero  when  dj^/dt  =  51/2  (for  x  ^0  andJTu 
7^0),  which  Implies  ^  changing  with  time.  There¬ 
fore,  the  separatrlx  arc  defined  by  equation  (l6) 
Is  no  longer  a  physically  possible  continuous 
trajectory  as  In  the  case  of  jl=  0,  but  only 
points  along  It  are  possible  Instantaneously  when 
d)^/dt  passes  through  fV2. 

Determination  of  the  Type  of  Singular  Points 

The  application  of  the ' techniques  described 
by  Coddlngton  and  Levinson  (2.)  to  equations  (8) 
reveals  the  following; 

1  The  critical  points  at  (+  1,0)  are  of 
higher  order  as  was  true  In  the  case  A  =  o,  but 
after  dividing  by  y  the  perturbations  f^^,  fall 
to  approach  zero  In  the  required  manner.  There¬ 
fore,  the  existing  theory  cannot  explain  their 
type.  However,  it  Is  known  from  equation  (6) 
that  these  €  phase-plane  critical  points  define 


the  half-frequency  whirl  phenomenon. 

2  The  singularities  at  (0,  +  1)  are  the 
"saddle  like"  type  Identical  to  those  at  the  same 
locations  for  the  case  of  A  =  o, 

3  The  new  critical  point  (also  the  equilib¬ 
rium  point)  at  X*  =  2/j  n  Is  a  center  for  the 
linear  case,  and  Is  either  a  spiral  point  or  cen¬ 
ter  for  the  nonlinear  case  (reference  7»  P*  382, 
Theorem  4.1).  Since  equation  (9)  is  symmetrical 
about  the  x-axls,  this  critical  point  must  be  a 
center  for  the  nonlinear  case  also.  A  center  Is 
symmetrical  —  a  spiral  point  Is  not.  Closed 
trajectories  indicating  neutral  stability  will 
surround  this  center  In  the  phase  plane.  Thus, 
if  the  journal  were  disturbed  slightly,  fhe  re¬ 
sulting  motion  would  be  continuously  oscillatory, 
corresponding  to  the  case  In  a  linear  system  of 
conjugate  complex  roots  on  the  Imaginary  axis . 

Journal  Behavior  for  Various  Values  of  J  A 

If  equation  (19)  Is  divided  by  JA,  and  jjv 
Is  considered  large,  the  conclusion  Is  that  x 
must  be  small  In  order  to  satisfy  the  resulting 
equation,  provided  (2/ri)  (d^/dt)  may  be  neglected 
In  comparison  with  one.  For  small  values  of  x, 
equation  (9)  becomes 

^  _  -  2  -  JOx  Vl  -  yZ  z  (24) 

dx  2  X  y 

Equation  (24)  has  singular  points  at  (0,  +  1)  and 
(2/jllO),  which  agrees  with  the  results  found 
previously  near  x  =  0.  Moreover,  the  singularity 
at  X  =  2/jAls  a  center,  and  trajectories  around 
the  center  reach  their  absolute  y  value  when  dy/ 
dx  =  0  in  equation  (24).  This  occurs  at  y  =  +  1. 


^  If  A  were  negative,  only  the  negative  x*  would 
satisfy  equation  (20),  but  only  positlveil -values 
are  considered  In  this  analysis.  If  the  effec¬ 
tive  rotational  speed  of  the  journal  were  In  the 
opposite  direction,  the  results  would  be  symmet¬ 
rical. 
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(a)  niase  Plane  (b)  Joiimal  Motion 

Fig.  7  Configurations  for  small  values  of  J  A 


Near  y  =  +  1,  pf  is  near  0  and  rr,  respectively. 

For  high  values  of  JjLthe  Journal,  If  disturbed, 
vfhlrls  at  very  small  eccentricity  ratios  with  at¬ 
titude  angles  continuously  changing  and  approach¬ 
ing  very  close  to  zero  and  -rr. 

Fig. 6(a)  shows  a  sketch  of  the  phase-plane 
portrait  taken  from  IBM  650  computer  solutions 
of  equation  (9).  The  left  half  plane  of  Fig. 6(a) 
Is  similar  to  that  for  the  case  of  fl.  =  0.  As 
trajectory  A,B,C,D,E,F  Is  traversed  In  Pig. 6(a), 
the  corresponding  physical  Journal  motion  Is  pic¬ 
tured  In  Pig. 6(b).  When  the  Journal  Is  located 
at  either  B  or  E,  any  slight  disturbance  due  to 
forces  not  considered  In  this  analysis  could 
cause  the  Journal  to  Jump  into  a  trajectory  lead¬ 
ing  toward  the  half -frequency  whirl  singularity 
at  a. 

Possible  confusion  as  to  the  Interaction  of 
the  negative  and  positive  half -phase  planes  In 
Fig. 6(a)  may  be  clarified  by  visualizing  the  neg¬ 
ative  half -plane  folded  about  the  vertical  y-axls 
onto  the  positive  half-plane.  This  type  of  men¬ 
tal  Image  Is  important  when  considering  cases 
where  trajectories  pass  close  to  (+1,  0)  due  to, 
for  example,  low  values  of  JA.  Pig. 7  Is  similar 
in  Interpretation  to  Pig. 6  with  the  exception 
that  the  possibility  of  travel  to  the  unstable 
singularity  occurs  at  an  eccentricity  ratio  near 
1. 

With  the  aid  of  equation  (3),  It  Is  seen  that 
large  values  of  JJl.  may  occur  because  of  high  ef¬ 
fective  Journal  speeds  and/or  light  Journal  loads. 
Other  factors  Influence  J,  such  as  viscosity, 
c/r,  and  so  on.  It  Is  also  worth  while  to  note 
that  since  the  Journal  attitude  angle  must  reach 
0  or  TT  for  the  Journal  to  pass  Into  the  unstable 
left  half -phase  region,  a  Journal  that  possesses 
both  large  J  (light  loading)  and  high  A- is  In  a 
particularly  vulnerable  position.  This  Is  be¬ 
cause  0  =  0,  V  at  y  =  J(dfe/dt)  *  ±  1,  and  for 
high  J  the  eccentricity  velocity  necessary  for 
o  0,  TT  Is  small.  This  deduction  gives  Impor¬ 
tant  Insight  Into  Journal  bearing  stability. 


SUMMARY  AND  CONCLUSIONS 

On  the  basis  of  the  assumptions  made  In  this 
analysis,  it  appears  possible  to  conclude  the 
following : 

1  Complete  Journal  bearings  with  continuous 
films  subjected  only  to  external  load  and  fluid 
pressure  forces  are  at  best  neutrally  stable. 

2  For  Journal  bearings  operating  at  high  ef¬ 
fective  rotational  speeds  or  large  values  of  J 
(light  loadings),  the  likelihood  of  unstable  mo¬ 
tion  appears  strong  If  the  Journal  Is  acted  upon 
by  a  sudden  disturbance.  This  Is  especially  true 
when  the  combination  of  both  high  speed  and  large 
J-values  exist. 

3  Unstable  motion  of  the  Journal  bearing  can 
also  occur  at  very  low  values  of  Jll. 
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APPENDIX  1 


Equations  (1)  and  (2)  are  of  the  form: 
cos 

sin  ^  =  g  -  k  ^ 


(25) 

(26) 


where  f,  g,  and  k  are  functions  of  t  and  con¬ 
stants.  Before  combining  equations  (25)  and  (26) 
and  eliminating  It  Is  convenient  to  restrict 
0  to  the  Interval  0  i  ^  and  allow  t  to  take 
on  any  of  the  values  In  the  Interval  -  1  ^  x  ^1, 
To  show  that  equations  (1)  and  (2)  are  unchanged 
by  this  restriction,  define 


0  when  0  4  4  * 

0  -  It  when  w  4  4  2  * 


,  '  _  *  when  0  4^4* 

-«  when  «  4  4  2  w 

where  0  -£  0  Zn  and  ^1.  If  0  0  -^-n , 

equations  (1)  and  (2)  are  clearly  unchanged.  If 
IT  ^  0  ^2ir.  then  cos  0'=  -  cos  0,  sin  0'=  -  sin 
0,  d^Vdt  =  -  df/dt,  d^Vdt  =  d^/dt,  and  0'=  -i  . 
Substitution  of  these  Into  equations  (1)  and  (2) 
yields  equations  In  and  t'  which  are  of  the 
same  form.  The  advantage  of  this  transformation 
Is  that  sin  0'  =  +  {1  -  cos^  0')^'^  and  the  am¬ 
biguous  minus  sign  Is  eliminated. 

Prom  equation  (25)  sin  0  Is  derived  from  the 


right  triangle  which  gives 

sjnr(  =  [1  -  (f  (27) 


If  equation  (25)  Is  differentiated  with  re¬ 
spect  to  time,  there  results 


-  sir  0 


A- 

d 


d  t 
d  t 


(28) 


Multiplying  both  sides  of  equation  (26)  by  sin  0 
gives 

sln^  0  =  g  sin  ^  -  k  sin  ^  (29) 


Substitution  of  equation  (27)  and  (28)  Into  equa¬ 
tion  (29)  and  rearranging  terms  gives 


k  { 


d2. 


d  t 


+  k 


d  f 
d  t 


d  « 

d  t 

+  g 


+ 

d  t 

\/l  - 
''  dt 


(50) 


-  1 


After  the  Indicated  differentiation  and  substitu¬ 
tions  are  performed,  equation  (50)  reduces  to 
equation  (7). 


APPENDIX  2 


Equation  ( 13 )  Is 


d  2  M  1  'd  t  ' 

(1  -  .V  (24  t 
2  € 


Let 

Then 


j  2  *  2 

( iJL\  =  «  =  V , 

'  d  t  '  , 

->  *2 
d^«  L**'  _  1  d  V 

7^2  "  2  d7-  -2d. 


(15) 


=  0 


Equation  (13)  becomes 

d  V  +  (7.^+2)  _  (  1  -  •  ^)^  (  2  +  ‘  ^  (31) 

d  .  .  (  1 


which  Is  a  linear  differential  equation  In  v 
whose  solution,  placea  In  pliase-plane  notation, 
reduces  to  equation  (15). 
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